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Abstract 

The  basic  task  of  a  battery  management  system  (BMS)  is  the  optimal  utilization  of  the  stored  energy  and  minimization  of  degradation  effects.  It 
is  critical  for  a  BMS  that  the  state-of-charge  (SoC)  be  accurately  determined.  Open-circuit  voltage  (OCV)  is  directly  related  to  the  state-of-charge 
of  the  battery,  accurate  estimation  of  the  OCV  leads  to  an  accurate  estimate  of  the  SoC.  In  this  paper  we  describe  a  statistical  method  to  predict 
the  open-circuit  voltage  on  the  basis  of  voltage  curves  obtained  by  charging  batteries  with  different  currents.  We  employ  a  dimension  reduction 
method  (Karhunen-Loeve  expansion)  and  linear  regression.  Results  of  our  modelling  approach  are  independently  validated  in  a  specially  designed 
experiment. 

©  2005  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

The  last  decade  has  seen  an  enormous  growth  in  number  of 
portable  devices  such  as  laptops  and  mobile  phones.  Develop¬ 
ment  of  such  devices  requires  a  drastic  improvement  in  power 
supply  and  power  management,  and  in  particular,  improvement 
in  the  capacity  of  rechargeable  batteries.  While  the  most  signifi¬ 
cant  progress  in  portable  energy  storage  depends  on  the  develop¬ 
ments  in  chemical  technology  and  material  sciences,  other  areas 
such  as  charge  regulation  and  battery  monitoring  can  contribute 
significantly  to  the  life-time  performance  of  these  batteries. 

It  is  of  most  importance  that  the  battery  management  system, 
which  controls  charging  and  discharging  of  the  battery,  operates 
with  an  accurate  estimate  of  the  energy  stored  in  the  battery 
at  any  given  time.  The  available  fraction  of  the  full  capacity  is 
called  the  state-of-charge  (SoC).  An  accurate  estimation  of  the 
state-of-charge  is  one  of  the  main  tasks  of  the  BMS,  Bergveld 
etal.  [1,2], 


*  Corresponding  author.  Tel.:  +31  620374232;  fax:  +31  402478190. 
E-mail  address:  snihir@eurandom.tue.nl  (I.  Snihir). 

0378-7753/$  -  see  front  matter  ©  2005  Elsevier  B.V.  All  rights  reserved, 
doi:  10. 1016/j  .jpowsour.2005 . 1 1 .090 


Several  conventional  approaches  exist  for  determining  the 
state-of-charge.  For  nickel-metal  hydride  (NiMH)  batteries,  one 
of  these  approaches  is  based  on  the  measurements  or  estimation 
of  electrical  properties  of  the  battery,  such  as  its  open-circuit 
voltage  (OCV).  By  definition,  the  open-circuit  voltage  is  the 
battery  voltage  under  the  equilibrium  conditions,  i.e.  the  voltage 
when  no  current  is  flowing  in  or  out  of  the  battery,  and,  hence 
no  reactions  occur  inside  the  battery. 

The  open-circuit  voltage  is  a  function  of  state-of-charge, 
OCV  =  /(SoC),  and  the  function  /  is  expected  to  remain  the 
same  during  the  life-time  of  the  battery,  i.e.  it  does  not  depend  on 
the  age  of  the  battery.  Note,  however,  that  other  battery  character¬ 
istics  do  change  with  time,  e.g.  capacity  is  gradually  decreasing 
as  a  function  of  the  number  of  charge-discharge  cycles.  In  this 
paper  we  only  consider  the  problem  of  estimating  the  OCV.  The 
question  of  how  exactly  the  OCV  depends  on  the  SoC  will  be 
addressed  elsewhere. 

In  this  paper,  we  study  the  voltage  of  the  battery,  which  is 
being  charged  with  a  constant  current.  Roughly  speaking,  with 
twice  the  current  the  battery  is  charged  in  half  the  time.  There¬ 
fore,  to  compare  the  results  of  charging  experiments  with  dif¬ 
ferent  currents,  we  study  the  battery  voltage  not  as  a  function  of 


1.  Snihir  et  al.  /  Journal  of  Power  Sources  159  (2006)  1484-1487 


1485 


time,  but  of  the  so-called  depth-of-charge  (DoC): 

It 

DoC(f)  =  — - , 

nom 

where  I  is  the  charging  current,  t  the  charging  time  and  (2nom  is 
the  nominal  capacity.  If  one  charges  an  empty  battery,  then,  typ¬ 
ically,  in  the  first  half  of  the  charging  period  the  state-of-charge 
and  the  depth-of-charge  coincide.  Later,  due  to  the  presence  of 
so-called  side-reactions,  the  SoC  increases  slower  than  the  DoC. 
Therefore,  in  the  figures  below  the  DoC  can  reach  values  larger 
than  100%,  while  SoC,  by  definition,  does  not  exceed  100%. 

The  open-circuit  voltage  corresponds  to  equilibrium  condi¬ 
tions.  Therefore  we  estimate  OCV  as  a  limit  of  voltage  curves 
corresponding  to  the  decreasing  sequence  of  currents.  More 
specifically,  using  voltage  curves  obtained  when  various  currents 
flow  into  the  battery  while  it  is  being  charged,  we  extrapolate 
to  the  current  value  equal  to  0,  hence  obtaining  the  open-circuit 
voltage  of  the  battery. 

The  second  approach  is  based  on  a  different  limiting  behav¬ 
ior.  If  a  battery  is  left  to  “rest”,  then  neglecting  some  physical 
processes  such  as  self-discharge,  we  can  assume  that  the  battery 
voltage  will  stabilize  to  a  certain  value,  which  is  the  value  of  the 
open-circuit  voltage  at  a  given  depth-of-charge.  The  relaxation 
period  can  take  some  time  and  the  OCV  is  determined  at  the  end 
of  the  resting  period. 

The  first  approach  uses  statistical  methods  to  predict  the  OCV. 
The  method  maps  a  voltage  curve  to  a  small  number  of  its 
Karhunen-Loeve  (KL)  coefficients;  each  curve  becomes  a  point 
in  a  low-dimensional  space  and  each  point  describes  a  voltage 
curve  obtained  for  a  specific  current.  A  polynomial  regression  is 
then  applied  to  model  the  dependence  of  the  KL  coefficient  on 
the  current.  Extrapolating  the  regression  model  to  zero  current 
yields  an  estimate  of  the  KL  coefficient  of  the  OCV  curve. 

The  second  approach  is  purely  experimental.  It  requires  a 
large  number  of  partial  charging  experiments  followed  with  sub¬ 
stantial  resting  periods.  This  approach,  however,  is  not  very  suit¬ 
able  for  the  application  in  battery  management  systems. 

A  good  agreement  between  the  results  in  OCV  estimation  by 
both  methods  is  demonstrated. 

All  the  experiments  were  performed  with  commercial  AA 
size  NiMH  batteries  of  the  type  HHR110AAOK  (Matsushita 
Battery  Industries,  Japan).  The  nominal  battery  capacity  was 
about  1100  mAh  at  room  temperature.  The  battery  was  placed 
on  a  special  holder  to  ensure  electrical  connections.  A  pre¬ 
calibrated  thermocouple  is  pasted  on  the  battery  to  monitor  the 
battery  temperature.  The  set  up  was  then  placed  in  an  ambient 
temperature  controlled  chamber.  The  battery  voltage  and  tem¬ 
perature  were  measured  during  the  experiment.  Electrochemical 
measurements  were  performed  using  Maccor  Systems  (Maccor 
Inc.,  Tulsa,  OK,  USA,  Series  2300).  Prior  to  all  experiments,  the 
battery  was  activated  according  to  a  standard  activation  proce¬ 
dure  to  ensure  battery  stability  and  hence  increase  reproducibil¬ 
ity  of  the  experiments  and  similarity  between  different  used  bat¬ 
teries.  The  activation  procedure  consisted  of  five  cycles  each 
containing  the  following  steps:  charging  with  —0.55  A  till  125% 
DoC  followed  by  a  discharge  phase  containing  a  rest  period  of 
1  h  then  discharging  with  —0.55  A,  a  rest  period  of  15  min,  fur¬ 


ther  discharging  with  —0. 1 1  A  and  a  final  rest  period  of  15  min. 
The  discharge  steps  have  a  battery  voltage  cut-off  value  of  about 
0.9V. 

2.  Open-circuit  voltage  estimation 

2.1.  Extrapolated  OCV  curve 

Let  us  consider  a  group  of  voltage  curves  V(t)  obtained  in  an 
experiment  when  a  NiMH  battery  was  charged  till  120%  DoC 
with  a  given  current  value  then  completely  discharged  according 
to  the  discharge  phase  described  in  the  activation  regime  and  then 
charged  again  with  a  different  current.  The  charging-discharging 
cycle  was  repeated  for  increasingly  charging  currents  ranging 
from  0.011  to  2.2  A. 

Since  the  measurements  were  performed  across  many  param¬ 
eters,  such  as  varying  currents,  voltages  and  time  points,  it  is  not 
possible  to  make  a  visual  inspection  of  the  relationship  between 
these  parameters  in  such  a  multi-dimensional  matrix.  One  way 
to  make  sense  of  this  data  is  to  reduce  its  dimension.  Several 
data  decomposition  techniques  are  available  for  this  purpose: 
Karhunen-Loeve  expansion  (a  form  of  principal  components 
analysis,  Flury  [3]  and  Jackson  [4])  is  one  of  such  techniques. 

Indeed,  the  representation  given  by  the  KL  expansion  with  a 
limited  number  of  KL  coefficients  is  optimal  in  the  mean-square 
sense.  Such  a  dimensionality  reduction  has  important  benefits. 
First,  noise  is  reduced,  as  the  data  not  contained  in  the  first 
coefficients  is  most  probably  due  to  noise.  Secondly,  in  some 
cases,  visualization  becomes  possible. 

Let  Lbe  an  n  x  m  matrix  representing  the  battery  voltage  dur¬ 
ing  charging  with  n  different  currents  and  sampled  in  m  points  in 
time.  We  assume  that  n  <  m .  Fig.  1  shows  a  typical  experimental 
result. 

The  singular  value  decomposition  (S  VD)  of  matrix  Y  is  given 
by 

Y  =  VZUT, 


Fig.  1 .  Voltage  during  charging  with  various  currents:  0.0 1 1 . 0.022, 0.055, 0.077, 
0.1 1,  0.55,  1.1,  1.65,  and  2.2  A  (bottom-up). 
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where  V  is  an  n  x  n  unitary  matrix,  S  —  diag(cri,  . . . ,  er„)  is 
an  n  x  n  diagonal  matrix,  with  cri  >  •  •  •  >  cr„  >  0,  and  U  is  an 
m  x  n  matrix,  whose  columns  are  the  first  n  columns  of  a  unitary 
matrix.  For  p  e  [1,  n],  the  Karhunen-Loeve  rank  p  approxima¬ 
tion  Yp  of  Y  is 

Yp  =  VpEpUj , 

where  Up  =  diag(<ri ,  . . . ,  ap ),  Vp  is  an  n  x  p  matrix  (first  p 
columns  of  V ),  and  Up  is  an  m  x  p  matrix  (hrst  p  columns  of 

U)- 

Karhunen-Loeve  approximation  Yp  of  Y  can  also  be  inter¬ 
preted  as  follows.  Denote  by  y [,...,  y„  e  R'"  the  rows  of  Y. 
Suppose  that  we  are  looking  for  a  linear  subspace  W  C  Rm , 
dim(lT)  =  p,  such  that 

n 

y II y*  -  Pw(yk)f 

k=  1 

is  minimal,  where  Pwiy)  denotes  the  projection  of  vector  y  to 
W.  Then  the  rows  of  L/J  form  the  basis  of  the  optimal  subspace 
W,  and  the  optimal  projections  are  given  by 

Pw(yk)  =  [oik,  1,  •  ■  • ,  oik,p]Uj,  k  =  1 - -  n. 

The  matrix  of  Karhunen-Loeve  coefficients  A  =  (o/k.  j)"/  ;=\ 
satisfies  A  =  VpSp. 

In  a  practical  application  of  the  KL  dimension  reduction 
method,  the  choice  of  the  parameter  p  depends  on  the  impor¬ 
tance  of  KL  basis  vectors  for  the  description  of  the  data,  which 
is  measured  by  the  associated  singular  values,  Karlis  et  al.  [5]. 
Typically,  for  the  data  sets  with  some  intrinsic  smoothness, 
most  of  the  original  variability  is  captured  by  the  hrst  few  basis 
functions.  For  our  experimental  data  the  hrst  three  Karhunen- 
Loeve  coefficients  seem  to  be  sufficient:  the  heuristic  indicators 
°aX.2XZ  =  °"  and  a/++-a+an  =  °-83  are  appropriately  large. 
Therefore,  we  put  p  —  3,  see  also  Fig.  2,  showing  the  hrst  three 
components. 

If  one  wants  to  construct  “experimental”  voltage  curves  that 
have  never  been  measured,  under  conditions  that  have  never 
been  used,  one  is  forced  to  interpolate  or  extrapolate  in  the 


DoC  [%] 

Fig.  2.  First  three  Karhunen-Loeve  basis  vectors. 


Current  [A] 

Fig.  3.  Dependence  of  Karhunen-Loeve  coefficients  on  current  value  in  low- 
and  high-current  regimes:  dotted  lines  represent  linear  regression. 

space  of  voltage  curves,  or,  equivalently,  in  the  space  of  KL 
coefficients.  Therefore,  one  has  to  model  the  dependence  of  the 
KL  coefficients  on  the  parameters,  influencing  the  experiment, 
e.g.  current,  temperature,  age  of  the  battery,  etc.  In  our  experi¬ 
ments,  the  value  of  the  applied  current  is  the  only  parameter  that 
has  been  varied.  Let  («/-  i ,  2,  oik. 3 )  be  the  KL  coefficients  de¬ 

scribing  the  voltage  curve  obtained  in  the  kth  experiment  when 
the  current  4  was  applied.  Fig.  3  shows  the  dependence  of  the 
first  KL  coefficient  041  on  4.  Fig.  3  is  quite  interesting  from 
the  electrochemical  point  of  view  as  well.  There  are  clearly 
two  different  slopes  in  the  dependence  of  oik ,1  on  h  f°r  l°w 
and  high  values  of  4.  respectively.  The  presence  of  two  differ¬ 
ent  slopes  is  related  to  the  fact  that  the  reaction  mechanism  of 
oxygen  evolution  has  two  rate-determining  steps,  see  Notten  et 
al.  [6]. 

In  this  paper  we  are  primarily  interested  in  the  open-circuit 
voltage,  i.e.  zero  current.  Therefore  we  apply  a  linear  regression 
model  for  the  KL  coefficients  using  only  low  currents  0.01 1— 
0.1 1  A,  i.e.  the  first  five  points  in  Fig.  3.  If  one  is  interested  in 
voltage  curves  corresponding  to  moderate  currents  (0. 1-0.5  A), 
kinetic  equations  derived  in  Notten  et  al.  [6]  provide  a  model  for 
interpolation. 

Since  we  adopted  a  linear  model  for  the  dependence  of 
(oik.i,  oik, 2,  oik, 3)  on  4  in  low-current  regime,  we  have  to  find 
only  six  regression  parameters  (ai,  b  1),  (122,  ^2)  and  (03,  £3), 
which  are  determined  by 

K 

T,  Wk.j  -  (a j  +  bjh) I2  min,  j  =  1,  2,  3, 

k=  1 

where  K  =  5  is  the  number  of  low  currents.  Finally,  setting  cur¬ 
rent  to  0  in  the  linear  model,  we  obtain  extrapolated  values  of  the 
KL  coefficients  a\  =  a\,  &2  =  (*3  =03,  and  the  prediction 

for  the  OCV  curve: 

A  /V  /V  /V  T 

V  —  [a ^2613]^  . 
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Fig.  4.  Partial  charging  experiment.  Battery  relaxation  period  is  6  h,  after  which 
the  OCV  is  sampled  (dots).  To  ensure  that  the  battery  reaches  equilibrium  state 
in  6  h,  a  very  low  current  has  been  used  for  charging. 


Fig.  5.  Open-circuit  voltage  as  a  function  of  DoC:  extrapolated  curve  (solid  line) 
and  several  experimental  values  (dots). 

2.2.  Experimental  OCV  curve 

The  values  of  the  OCV  for  various  values  of  the  depth-of- 
charge  can  be  measured  in  a  separate  experiment.  Comparing 
these  results  with  the  results  of  our  statistical  approach  we  can 
independently  validate  our  modelling  assumptions. 

In  this  experiment,  the  battery  is  charged  with  a  constant  cur¬ 
rent  to  a  specific  depth-of-charge,  the  current  is  then  interrupted 


and  the  battery  is  allowed  to  rest  for  a  certain  period  of  time.  The 
basic  idea  is  that,  after  current  interruption,  the  battery  reaches 
the  equilibrium  state  (6h  in  our  experiment)  and  hence  the  bat¬ 
tery  voltage  is  the  value  of  the  OCV  at  the  given  DoC.  Fig.  4 
shows  the  voltage  curves  obtained  in  such  experiments  for  five 
values  of  DoC.  Fig.  5  shows  a  good  agreement  between  the  re¬ 
sults  of  our  data  analysis  (solid  line)  and  several  experimentally 
observed  values  of  OCV  (dots). 

The  experimental  approach  to  the  evaluation  of  the  OCV  is 
extremely  time-consuming.  Every  point  on  the  OCV  curve  is 
obtained  by  performing  a  separate  charging  experiment.  On  the 
other  hand,  the  modelling  approach  allows  a  simultaneous  esti¬ 
mation  of  the  OCV  for  all  values  of  the  DoC  based  only  on  a  few 
charging  experiments  with  various  currents.  This  aspect  of  our 
approach  makes  its  application  in  adaptive  battery  management 
systems  feasible. 

3.  Conclusion 

In  this  paper  we  propose  a  method  for  determining  the  open- 
circuit  voltage.  Our  approach  is  based  on  Karhunen-Loeve  di¬ 
mension  reduction  and  modelling  the  dependence  of  the  KL  co¬ 
efficients  on  physical  parameters  of  the  experiment.  The  method 
is  also  suitable  for  interpolation  in  the  space  of  experimental 
voltage  curves.  Strikingly,  results  of  our  statistical  analysis  also 
detect  an  interesting  electrochemical  phenomenon  recently  ob¬ 
served  for  Ni-based  rechargeable  batteries.  We  also  performed  a 
series  of  independent  experiments  in  which  the  OCV  was  eval¬ 
uated  directly  for  several  values  of  the  DoC.  Results  of  both 
approaches  are  in  a  good  agreement. 

The  method  presented  in  this  paper  suggests  a  significant  po¬ 
tential  of  the  dimension  reduction  techniques  in  the  development 
of  battery  management  systems. 

References 

[1]  HJ.  Bergveld,  W.S.  Kruijt,  P.H.L.  Notten,  Electronic  network  modelling  of 
rechargeable  NiCd  cells  and  its  application  to  the  design  of  battery  manage¬ 
ment  systems,  J.  Power  Sources  77  (1999)  143-158. 

[2]  HJ.  Bergveld,  W.S.  Kruijt,  P.H.L.  Notten,  Battery  Management  Systems: 
Design  by  Modelling,  Kluwer  Academic  Publishers,  Boston,  2002. 

[3]  B.  Flury,  Common  Principal  Components  and  Related  Multivariate  Models, 
Wiley,  Chichester,  1988. 

[4]  J.E  Jackson,  A  User’s  Guide  to  Principal  Components,  Wiley,  New  York, 
1991. 

[5]  D.  Karlis,  G.  Saporta,  A.  Spinakis,  A  simple  rule  for  the  selection  of  principal 
components,  Commun.  Statist.  Theory  Meth.  32  (2003)  643-666. 

[6]  P.H.L.  Notten,  E.A.  Verbitskiy,  W.S.  Kruijt,  HJ.  Bergveld,  The  oxygen  evo¬ 
lution  and  recombination  kinetics  in  sealed  rechargeable,  Ni-based,  batter¬ 
ies,  J.  Electrochem.  Soc.  152  (2005)  A1423-A1433. 


